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' This paper completes investigation of symmetry properties of nonlinear variable coefficient 

, diffusion-convection equations of the form f(x)ut = {g{x)A(u)ux)x + h{x)B{u)ux which was 

^SJ ' started in [9-lf]. Potential symmetries of equations from the considered class are found and 

the connection of them with Lie symmetries of diffusion-type equations is shown. Exact 
(— I ' solutions of the Fujita-Storm equation ut = {u~'^Ux)x are constructed. 

3 ' 1 Introduction 

This paper is the fourth part of a series of works on symmetry properties of nonlinear variable 
coefficient diffusion-convection equations of the form 

^ ■ f{x)ut = {g{x)A{u)ux)x + h{x)B{u)ux, (1) 

'sj- ■ where / = f{x), g = g{x), h = h[x), A = A[u) and B = B{u) are arbitrary smooth functions of 

Q . their variables, f{x)g{x)A{u)^Q. 

In the first three parts of the series we studied equivalence transformations of equations ([T]), 
performed the complete group classification of class ([T]), completely classified local and poten- 
tial conservation laws of the considered equations, found contractions of their Lie symmetry 
algebras, solutions, conservation laws and characteristics of conservation laws. Using different 
^ ' generalizations of equivalence of conservation laws proposed in [11], we constructed a complete 

■ list of inequivalent potential systems of equations ([T]). Now, in the last part of the series we clas- 

sify potential symmetries of equations ([1]), i.e., symmetries of the constructed potential systems 
that are not projectible on the space of non-potential variables {t,x,u). 

The notion of potential symmetries was introduced in [4,5]. A system of partial differential 
equations may admit symmetries of such sort when some of its equations can be written in 
a conserved form. After introducing potentials for equations written in the conserved form as 
additional dependent variables, we obtain a new (potential) system of differential equations. Any 
local invariance transformation of the obtained system induces a symmetry of the initial system. 
If transformations of some of the "non-potential" variables explicitly depend on potentials, this 
symmetry is a non-local (potential) symmetry of the initial system. 

Potential symmetries of different subclasses of class ([1]) have been investigated by many 
authors. Thus, e.g., in [4] potential symmetries of constant coefficient diffusion equations 
{f = g = 1, B = 0) were found, potential symmetries of equations with f = g = h = 1 
were found in [26]. In [23] potential symmetries of constant coefficient equations ([T]) were clas- 
sified with respect to a group of potential equivalence transformations and connection between 
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potential and Lie symmetries was found. In [17] a number of results concerning (local and poten- 
tial) equivalence transformations of constant coefficient equations ([T]) were obtained. Potential 
symmetries of variable coefficient diffusion equations {B = 0) were classified in [27-29]. 

Note 1. Since potential symmetries of linear parabolic equations have been studied in [25], in 
the subsequent analysis we consider only nonlinear equations of form ([1]). 

The paper is organized as follows. First, we briefly recall results on equivalence transfor- 
mations and a list of inequivalent potential systems (Section [2]) of equations from class ([T]). 
Preliminary analysis of the derived potential systems in quite a general form is given in Sec- 
tion [3l The complete classification of potential symmetries of equations from class ([T]) is given 
in Sectional where we also describe a connection between potential and local symmetries of ([T]). 
Although Fujita-Storm equation is a model for a wide range of physically important processes, 
only few of its exact solutions are known. Therefore, in Section ([5]), applying potential equiv- 
alence transformation, we adduce some new exact solutions of the Fujita-Storm equation that 
cannot be obtained by usual Lie symmetry reduction method. 



2 Equivalence transformations and potential systems 

The usual equivalence group G~ of class ([T]) is formed by the nondegenerate point transforma- 
tions in the space of {t, x, u, /, g, h, A, B), which are projectible on the space of (t, x, u), i.e. they 
have the form 

{i, X, u) = {T\ T^, T-){t, X, u), 

if, g, h, A, B) = {Tf, ^^ T^, r^)(t, X, u, f, g, h, A, B), (2) 

and transform any equation from class ([1]) for the function u = u{t, x) with the arbitrary elements 
(/, g, h, A, B) to an equation from the same class for function u = u{t, x) with the new arbitrary 
elements {f , g, h, A, B). 

Theorem 1. [8, 9] C consists of the transformations 
t = 5it + 52, x = X{x), u = 53U + 54, 

f = ^r'f, 9 = £i£2^X^g, h = ei£^^h, i = £2-4, B = e2,B, 

where 5j {j = 1,4) and £i {i = 1,3) are arbitrary constants, 5i6s£i£2£3 ^ 0, X is an arbitrary 
smooth function of x, 7^ 0. 

Class ([1]) admits other equivalence transformations which do not belong to G~ and form, 
together with usual equivalence transformations, an extended equivalence group. We demand for 
these transformations to be point with respect to (t, x,m). The explicit form of the new arbi- 
trary elements (/, g, h, A, B) is determined via (t, x, /, g, h, A, B) in some non-fixed (possibly, 
nonlocal) way. 

Theorem 2. [8,9] The complete extended equivalence group is formed by the transformations 
i = 5it + 82, X = X(x), u = 5^u + 54,, 

f = -i^f, g = eie^^X^ipg, h = eie^^(fh, i = £2^, = £3(5 + £4^), 
where Sj {j = 1,4) and ei {i = 1,4) are arbitrary constants, 5i(53£i£2£3 / 0, X is an arbitrary 

-SA f H^dx 

smooth function of x, X^ ^ 0, (f = e . 
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Group contains a (normal) subgroup of gauge equivalence transformations 

f = eiipf, g = £ie2^ipg, h = £ie^^(ph, A = £2A, B = £3(3 + £4A), (3) 

where = e ^'^■^ ai^)^^^ £■ [i = 1^4) are arbitrary constants, £i£2£^ 7^ (the variables t, x and 
u do not transform!). 

The transformations ([3]) act only on arbitrary elements and do not really change equations. 
Application of "gauge" equivalence transformations is equivalent to rewriting equations in an- 
other form. 

The factor-group /G^^ coincides for class ([T]) with G^ /G^^ and can be assumed to consist 
of the transformations 

t = 5it + 82-, X = X{x), u = 53n + (54, 

/ = ^/, ~g = X^g, h = h, A = A, B = B, 

where 5i {i = 1,4) are arbitrary constants, 8163 7^ 0, X is an arbitrary smooth function of x, 
X, / 0. 

Using the transformation i = t, x = f ^|y, u = u from G^/G'^^, we can reduce equation ([1]) 

to 

/(x)% = {A{u)us;)i + h{x)B{u)us:, 

where f{x) = g{x)f{x), g{x) = 1 and h{x) = h{x). (Likewise any equation of form ([T]) can be 
reduced to the same form with f[x) = 1.) That is why, without loss of generality, we can restrict 
ourselves to the investigation of the equation 

f{x)ut = {A{u)ux)^ + h{x)B{u)ux. (5) 

Any transformation from G"", which preserves the condition (7 = 1, has the form 

i=5it + 52, X = 65 f e^^ f'dx + 65, u = d3U + 64^, 

f = 6i6^^6gfe-^^^ I\ h = 6g6:f^he-^^ /^ (6) 
i = 55M, B = d-r{B + 68A), 

where 5i {i = 1,9) are arbitrary constants, SiSsS^djSg ^ 0. (Here and below Jh = Jh{x)dx.) 
Note 2. If i? = 0, we assume h = I. 

Note 3. It has been shown in [11, 24] that the equivalence group for a class of systems or 
the symmetry group for single system can be prolonged to potential variables for any step of 
the direct iteration procedure. It is natural to use the prolonged equivalence groups to classify 
possible conservation laws and potential systems in each iteration. Additional equivalences 
which exist in some subclasses of the class or arise after introducing potential variables can be 
used for further analysis of connections between conservation laws. 

In the third part [11] of the given series we found all possible inequivalent (with respect to G~ 
and symmetry transformations of the corresponding classes of equations) simplest potential 
systems of nonlinear equations 

1. h = l: Vx = fu, vt = Au^ + J B 

2. A = l, Bu / 0, / = -h{h-^)^,: = e*(/i-i),,^/, vt = e\-h-^u^ + {h-^)-,u - JB). 
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3. B = 1, f = h^: Vx = e^h^u, Vt = e^{Aux + hu). 

4. B = 1, f = hx + hx~^: Vx = e^xfu, vt = e^{xAux + xhu — JA). 

5. B = 0: Vx = xfu, Vf = xAux — J A 

6.1. B = 1, f = h = e-^/^x~\ fi G {0, 1}: 

Vx = e^^^xfu, Vt = e'^^x{Aux + hu) — e'^* jA, 

6.2. B = 1, f = h = e-''/^x-^ fi G {0, 1}: 

Vx = e''*(tx - l)/n, Vt = e^*(te - l){Aux + hu) - te^* jA 

7.1. B = l, f = \x- l|'^-3/2|x + Ij-A'-s/s, /j = 1^ _ i|/^-V2|^ + l|-/^-i/2. 

Va; = e2'^*(xcosht — smht)/^, ft = e'^'^^{xcosh.t — smh.t){Aux + /in) — e'^'^^ cosht jA. 

7.2. 5 = 1, / = |x - l|'^-3/2|x + ij-A'-s/a^ = |x - Ij/'-^/^i^ _^ l|-/^-V2. 

= e2'^*(xsmht — cosh.t)fu, vt = e^'^^{xsinht — cosht){Aux + hu) — e^'^^sinht jA. 

7.3. B = 1, f = \x- l|'^-3/2|x + h = \x- l\^-^/^\x + l|-/^-i/2. 

= e(2/^+i)t(x - l)fu, Vt = e(2/^+i)*(x - l){Aux + /i-u) - e^^A^+i)* jA, 

8. B=l, f = e/^'^'-'=t'^'i^(x2 + l)-3/2^ /i = gMarctanx(^^2 _^ l)-l/2. 

= e^*(xcost + smt)fu, vt = e^*(xcost + sint)(Ana; + /in) — e^* cost jA. 

For some subclasses of ^ the spaces of local conservation laws are two-dimensional, that 
allowed us to construct general potential systems by means of introducing two potentials for 
each case simultaneously: 

5'. B = Q: Vx = fu, Vt = Aux, Wx = xfu, wt = xAux — JA. 
6'. B = l, f = e-Z'/^x-^, h = e-^/^x-i, fi € {0, 1}: 

Vx = e^^xfu, Vt = e'^^x{Aux + hu) — e^* JA 

Wx = e^\tx - l)fu, Wt = e^\tx - l){Aux + hu) - te^* jA. 

7'. B=l, f = \x- l|/^-3/2|x + l|-A'-3/2^ h=\x- l|/^-^/2|x + l|-^-^/2: 

Vx = 6(2^^+1)* (x - l)fu, Vt = e(2/^+i)*(a; - l){Aux + hu) - e(2M+i)t jA 
Wx = e(2A'-i)*(x + l)/n, Wt = e(2A'-i)*(x + l){Aux + hu) - e(2M-i)* jA. 

8'. B =1, f = Cretan a' ^^2 _^ l)~3/2^ ^ g/. arctan x ^^2 ^ l)-l/2. 

Vx = e^*(xcost + smt)fu, vt = e^*(xcost + sint)(Ana; + hu) — e^* cost jA, 
Wx = e'^*(x sint — cos t)fu, wt = e^*(x sint — cos t)(Aux + hu) — e^* sint jA. 

For more details on equivalence relations on the set of potential systems and rigorous defini- 
tions and proofs see [11]. 

In the subsequent analysis we classify potential symmetries of equations ([5]) arising from the 
above potential symmetries with respect to (trivial) prolongation of transformations from 
to the corresponding potential variables. 
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3 Preliminary analysis of potential systems 



As one can see, all simplest potential systems of equations ([5]) have the similar structure: 

Va: = a{t,x)u, vt = P{t,x)A{u)ux +'jit,x,u), (7) 

where af3 ^ 0. Here a{t, x) is nothing but the characteristic of the corresponding conservation 
law. 

The initial equation on u is a differential consequence of system ([7]). Another differential 
consequence is the equation 

vt = Pr-) +7(i,^,-) (8) 

on the potential dependent variable v which is called the potential equation associated with the 
equation ([5]) and the characteristic a. 

Consider a Lie symmetry operator Q = rdt + ^dx + V^u + Gd^ of system d?]). The coefficients 
of Q are functions of t, x, u and v. The infinitesimal invariance criterion implies for system d?]) 
the following system of determining equations 

Tu = ^u = Ou = 0, (9) 

rx = r^ = 0, (10) 
r? = ^ l^-S.va'^u'^ - {atT + a^x + "xC - a6y)u + 9^^ , (11) 

= -Tlu + Ov+ix- — -Tt- — , (12) 

-PA{r]x + ari^u) + j{0y - n - a^vu) - 7«r? - 7xC - itT - a^tu + 9t = 0. (13) 

The subsystem ([9]) of determining equations means that any Lie symmetry transformation 
of d?]) with respect to t, x and v does not depend on u. Equation ([8]) is a differential conse- 
quence of system ([7|), and there exists one-to-one correspondence between the sets of solutions 
of equation ([8]) and system Therefore, the truncated operator Q = rdt + + Od^ is a Lie 
symmetry operator of equation ([8]). 

And vice versa, consider a Lie symmetry operator Q = rdt + ^dx + 99^ of equation ([8]). The 
coefficients of Q are functions of t, x and v. Then the prolonged to u operator Q = Q + r]du, 
where r] is defined by formula (jlip . is a Lie symmetry operator of system ([7]). 

General potential systems of equation ([5]) look like 

Vx = a{t, x)u, vt = (3{t, x)A{u)ux + 7(t, x, u), 

Wx = X{t,x)u, wt = fJ-{t, x)A{u)ux + I'it, x,u), (14) 

where a(3Xfj, ^ 0, a and A are linearly independent. 

Applying Lie infinitesimal criterion for system (jl4p being invariant with respect to the sym- 
metry generator Q = rdt + S,dx + r]du + 98^ + Cdw, similarly to the case of the simplest potential 
system, we get 

"^u — — 9u — Cu — 0, 

Tx — Ty — Tii) — 0. 

All potential systems of forms ([7]) are ()14p are constructed with usage of local conservation 
laws of equations ([5]) which were classified with respect to G^- Each of these subclasses of 
conservation laws is invariant with respect to some subgroup of ■ each case this subgroup 
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has a very simple structure and can be singled out from Q imposing additional condition 
^7 = ^8 = 0. It becomes the usual equivalence group of class ^ and is trivially prolonged to 
potentials. Henceforth we will call such prolongation as G^j.. 

More detailed analysis of determining equations for coefficients of symmetry generators for 
both simplest and general potential systems requires direct substitution of arbitrary elements a, 
f3, 7, /X, u and tp. We skip the intermediate cumbersome calculations and adduce immediately 
the final results. Namely, in the next section we list all possible G^-inequivalent pure potential 
symmetries of equations ([5]), i.e., symmetries of the derived potential systems that are not 
projectible on the space of non-potential variables {t,x,u). 



4 Potential symmetries and connections with local symmetries 

Now we will investigate Lie symmetries of each of these potential systems separately. Below 
we analyze these symmetries adducing only the cases when the Lie symmetries of the systems 
induce nonlocal symmetries of the initial equation. At first, let us state that potential systems 
2, 6.1-7.3, 6'-8' do not yield any potential symmetry for equations ([5]). 

Potential system 1 associated to the simplest characteristic 1 arises in case h = 1 and, as we 
pointed out above, has the form 

Vx = fu, vt = Aux + J B. (15) 

Its constant coefficient case (/ = 1) has been studied in [23,26]. In particular, ibid we performed 
the complete group classification of system (fT5]) |/-i and proved that all the potential symmetries 
can be obtained from Lie symmetries of equations (l5])|j=i by means of prolongation to the 
potential v and application of potential equivalence transformation 

_ u 

t = t, X = X + ev, 



1 + eu 

and potential hodograph transformation 

i = t, X = V, u = u^^, V = X. (16) 



That is why for the constant coefficient case of system (|15p below we adduce only the list of 
systems and their Lie symmetries which induce potential symmetries of equations from class ([5]) 
and do not consider their connection with Lie symmetries of diffusion equations. Moreover, 
we assume that connection of a potential symmetry of a variable coefficient equation to Lie 
symmetry is found by reducing the corresponding potential system to a system (fT5]) | j-i . 

The list of G^-inequivalent constant coefficient systems (jlSp inducing potential symmetries 
of equations from class ([5]) consists of the following ones: 

1. f = l,A = ^x-2eV«, B = 0: A^^^^ = {dt, dx, d,, 2tdt + xdx + vd,, tdt - vd^ + u^du); 
^max ^ (^Q^^ 2tdt + xdx + vd^, ^itdt + vdx - u{u + l)du - vdy); 

/i arctan u 

3. f = l,A= B = 0: 

+ I 

^max ^ (^Q^^ 2tdt + xdx + vdy, ^tdt + vdx - {y^ + 1)5„ - xd^); 

4. / = 1, ^ = S = 0: 

^max ^ i^Q^^ 2tdt + udu + vdv, -vxdx + u{ux + v)du + 2tdy, 

At^dt — {v^ + 2t)xdx + n(v^ + 6t + 2xuv)du + Atvd^, xdx — udu, 4>dx — (kvU^du)] 
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5. / = 1, A = li-SeA*/", 5 = - l)ei/": 

^max ^ (^g^^ g^^ g^^ _ 2)tdt + {{fl - l)x + v)d, - U^O^ + (/i - l)^;^^); 

6. f = 1, A = ti-^e^/", B = -u-^: 

^max ^ ^g^^ g^^ g^^ ^g^ + + _ ^2g^ + _ 2^)^^). 



^max ^ ^g^^ 9„ (;u - 2z.)tat + ((^ - - v)d^ + n(u + + l)vd,); 
8. / = l,^=_^i-_S = -(n+l)ln 



('U+l)/^+2 n+1 u + 

^max ^ ^0^^ ^^g^ + (^^ + ^ _ t)^^ + ^(^ + l)^^ + + l)^;^,); 

„/iarctanu _i_ , .\ „i/ arctan m 

9. / = i, ^ = 1^ _ ^x + ^J^ , 

^max ^ ^5^^ g^^ _ 2^,)^^^ + ((^ _ j,)^; _ + (^2 ^ + (a; + _ u)v)d^)- 

10. / = 1, ^ = n-^ 5 = -n^2. 

^max ^ l^Q^^ Q^^ 2tdt + + Vd^, -Vd^ + U^du + 2td^, 

W^dt - (f^ + 2t)d^ + 2u{uv + 2t)du + ^tvd^, dx, e-^(/)5^. + e-^(0 - 

11. f = 1, A = l, B = 2u: 

^max ^ ^^^^ 2tdt + - na„, 2tdx -du- xd,, 

W^dt + Atxdx - 2{x + 2ut)a„ - (x^ + 2t)5^,, 5^,, e'^iK - hu)du + e-^'hd^). 

Here = const. i^) 7^ (—2,-2), (0,1) and 1/ / —1,0 for case 7. /i 7^ —2,0 for case 2. 
The functions (j) = <j){t, v) and h = h{t, x) are arbitrary solutions of the linear heat equation 
{(j)t = (j)vv; ht = hxx)- 

Up to the group G~ there exists exactly one system of form (fT5]) with non-constant value 
of / yielding potential symmetries of equation ([5]), namely, 

Vx = X~^^^U^ Vt = U~'^Ux (17) 

associated to the diffusion equation x~^/^ut = {u^'^Ux)x- The corresponding algebra of potential 
symmetries is 

{dt, dy, 3xdx — udu — 2vdv, 2tdt + udu + vdy, 3xvdx — {v + 3x~^^^u)udu — v^d^). 

Application of potential hodograph transformation maps system ()17p to the system 

u = vi = {v~'^/'^Vx)x- 

Therefore, we establish the nonlocal mapping between the equations x~'^/'^ut = {u~'^Ux)x and 
vi = {v~^^^Vx)x belonging to class ([5]). (In implicit form this transformation was written firstly 
in [18].) The same transformation establishes isomorphism between the above algebra of the 
potential symmetries and maximal Lie invariance algebra 

{dti dx, 2tdt + xdx, 4:tdt + 3^9^,, x'^dx — Sxvdy) 

of the equation Vf = {v~^^^v^)x- Tildes over variables in the Lie symmetry operators are omitted. 

Considering this together with results of [23, 26] we can state that potential symmetries of 
equations dS]) obtained from system (|15p are completely classified and connection between the 
potential symmetries and Lie symmetries of equations ([5]) are found. 

Potential symmetries obtained from system 3. It is possible to introduce potential system 
associated with the characteristic e* of form 

Vx = e^hxU, Vt = e^{Aux + hu). 
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for subclass B = 1, f = of class ([5]). There exist five G^-inequivalent values of arbitrary 
elements for which system 3 gives potential symmetries of the initial equations. The first four 
of them are quite similar: 

1. / = 1, /i = X, ^ = n-^e"^/": 

2. f = 1, h = X, A = : 

^max ^ (e-2*(at-x5,), e'^S,, d,, dt + vd,, -dt + {fix-2e-'v)d, + 2u{u+l)du + {fi + l)vd,); 



e 



jj, arctan u 



3. f = l,h = x,A 

A^-''- = {e-^\dt-xd^), d,, dt + vd,, {fix - 2e-'v)d^ + 2{u^ + l)du - {fiv - 2e'x)d,); 

4. f = 1, h = X, A = u-^: 

^max _ (^^~2t^Q^ _ xdx), dy, xdx — udu, dt + vd^, —vxdx + {v + e^xu)udu + e^*9„, 
e^^Ot - (2e2* + v'^)xda: + (Se^* + 2e^xuv + v'^)udu + 2e'^^vdy, e'^ifd^ - (fvu'^du). 
Here ip = ^{t,v) is an arbitrary solution of equation e~'^^ipt = ipw 

Using the local transformation i = e'^^/2, x = xe*, u = u we can reduce them, and moreover, 
any equation of form Uf = {A{u)ux)x + xux to the constant coefficient equation -Uj = {A{u)ux)x- 
The same transformation trivially prolonged to the potential v = v maps the corresponding 
potential system 

Vx = e^u, vt = e^{Aux + hu) 

to the constant coefficient potential system 

Vx = U, Vi = A{u)ux 

associated to the characteristic 1 (potential system 1) that have been studied in [4,23,26]. 
The last nontrivial example of the systems 3 giving potential symmetries of equations ([S]) is 

5. f = hx, h = x~^/^, A = u-'^: 

^max ^ (^dt + vdy, dy, 2tdt + 3xdx-udu, e^\dt + 5xdx), 3xvdx + u{e^x-^/^u-v)du-v^d^). 
These values of arbitrary elements correspond to the equation 

-x~^/^nt/3 = {u~'^Ux)x + x-^^^Ux. 

Using the point transformation i = 3(e^^* — 1)/2, x = xe~^*, u = u one can reduce it to equation 

X'^^/^U^ = {u~'^Ux)x 

The same transformation trivially prolonged to the potential variable v = v maps the corre- 
sponding potential system to (fT7|) . considered earlier. 

Potential symmetries obtained from system 4 for subclass B = 1, f = hx + hx^^ of 
class ([5]) are ones of the most interesting in the sense of reducibility of them to the local sym- 
metries. Namely there exist two inequivalent systems of form 4 giving potential symmetries of 
equations 

1. h = x^/^, A = u-'^: 

^max ^ {dt+vdy, dy, 3xdx-2udu+2vdy, e-^{4dt-3xdx), 3xvdx-2u{v+2e^x^/^u)du+v'^dy). 
The point transformation i = (e^*/^ — l)/2, x = xe^*/^, u = u maps this equation to 

x~'^/^ui = {u~'^Ux)x 

and prolonged transformation v = v maps the corresponding potential system to system 5.1 
(ci = 0, C2 = (4/3)~^/^) that will be considered below. 
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2. h = \n \x\/x, A = u A"^^"^ = {dt + vdy, e ''xd^, e ^xvd^ — u^du + 2e*(9„) 
Transformation t = t,x = v,u=l/u, v = x maps the given potential system to the 

— = e M, — = e — In f — It , 

that is a potential system of generalized Burgers equation e~^ui = e^Uxx — "^uux- Transforming 
additionally time variable as t = e* we obtain generalized Burgers equation 

studied at [14] with three-dimensional symmetry algebra {dx-, idf. + xdx, idx + ^du). Note, that 
the projection of the operator di, on the space of variables {i,x,u) is trivial. 

Potential symmetries obtained from system 5. Consider the last simplest potential system 
Vx = xfu, vt = xAux — J A (18) 

associated to the characteristic x corresponding to the subclass = of equations ([5]). The 
following two cases correspond to the nontrivial potential symmetries of the initial equation. 

1. / = (cix3/2 + C2xV2)-4/3^ = ^-2. ^max ^ 

{dt, dy, 4c2tdt-3{cix + C2)xdx + {4:C2 + 3cix)udu, 3{cix + C2)xdx-{2c2 + 3cix)udu + 2c2vdjj, 
3{cix + C2)xvdx — {2>x^/^{cix + C2)^^^'^u + (2c2 + 2>cix)v)udu + C2v'^dy) 
Usual potential hodograph transformation (jl6p maps arbitrary system of form (jlSp to 

u 

Vx = -r-^, Vl = VUx-UVx- 

vf 

Eliminating variable u we get a constant coefficient equation —{\/v)i = {f{v)vx)x- After chang- 
ing the variable v = 1/v we obtain constant coefficient equation 

Vi = {f{v)Vx)x 

having the form ([5]). Here / = v^'^f{l/v). Application of this transformation to system with 

/ = (cix3/2 + C2a;l/2)-4/3 Igads to / = {C2V + Ci)-^'^ . 

2. f = x-'^, A = u-^: A"^'''' = {dt, dy, xdx,2tdt + udu + vdy, xvdx-u'^du + 2tdy, 

it'^dt + (u^ + 2t)xdx + (4t — 2uv)udu + Atvdy, x^kpdx — xu{ip + (pyu)du). 
Here if = (p{t,v) is an arbitrary solution of equation (ft = ifyy. 

Local transformation i = t, x = —1/x, u = xu, v = v maps the equation x~'^Ut = {u~'^Ux)x 
into Uf: = {u~'^Ux)x and the given potential system to 

Vx = u, vi = ir'^Ux, 

which is potential system of the constant coefficient equation associated to the characteristic 1. 
Potential symmetries obtained from system 5': i? = 0: 

Vx = fu, Vt = Aux, Wx = xfu, wt = xAux - JA. (19) 

Under the transformation v = v, w = xv — w this system is locally equivalent to the second- 
level potential system 

Vx = fu, Wx = v, wt = jA (20) 
that can be obtained from the conservation law Dt{v) — Dx{jA) = of the potential system ([15]) . 
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Below we adduce classification result for system (j20p omitting the tildes over the variables 

1. / = ^ = n-2/3: 

{dt, dy + xdji,, dyj, 2tdt — 3xdx — I5udu — 3ti;5^, Atdt — 3xdx — I2udu + 3vdy, 
x^dx + Sxudu + {w — xv)dv + xwdw, xwdx — ?>{xv — 2w)udu — v{xv — w)dv + w^dw)] 

2. / = 1,^ = ^-2/3; 

{dt, dx, dv + xdw, dw, 2tdt + xdx + vdy + 2wdw, 2tdt + 3udu + 3vdy + Swd^, 
wdx — Suvdu — v^dy); 

3. f = x-'^, A = u-^: 

{dt, dw, xdx - vdy, 2tdt + xdx + udu + wd^, 
x{2xv — w)dx — u{xv + 2u)du + v{w — xv)dy + (x^f^ — 2t)dw, 
At^dt + x(6t + 3x^1!^ + -Axvw + w'^)dx + 2u{2t - 3xuv + 2uw - x^u ^ + xvw)du 
+ {2xvw — 2t — x^u^ — w'^)dy + 2(2t?i; + x^v^ — x'^v'^w)dw, 
x^Pzdx + (Pz - upzz)xudu + pdy + x{p + xvpz)dw). 
Here z = if — xf , p = z) is an arbitrary solution of the linear heat equation pt = Pzz- 

4. / = 1,^ = ^-2; 

{dt, dv + xdw, 2tdt + xdx + vdy + 2wdw, xdx - udu + wdw, 

{w — 2vx)dx + {2xu^ + uv)du + 2tdu + (2tx — xv'^)dw, X^dx — \tv?du + {vXv — X)dw, 
2fdt + {vw - |xt;2 - 'itx)dx + {u^{'ixv -w) + ^ + 5tu)du + 2tvdy + {^v'^w -tw- xv^)du,) 
Here X = X{t , f) is an arbitrary solution of the equation — Xy^. 

5. / = 1, >1 = (u2 + l)'i: 

(Ot, dx, dy+xdw, dw, 2tdt + xdx + vdy + 2wdw, -vdx + {l + u'^)du + xdy + ^{2t + x'^ -v'^)dw); 

6. / = 1, ^ = (n2-l)-i: 

Si:, 9^^,, dy + xdw, 2tdt + xdx + vdv + 2wdw, vdx + xd^- {u'^ -l)du + {-t+ ^ + \)dw); 

7. / = 1, yl = (n2 + l)-ie^^^'=*^'"": 

dt, dx, dy + xdyj, dyj, 2tdt + xdx + vdy + 2wdu,, p.tdt-vdx + {l + v?)du + xdv + {x^ -v'^)du,); 

8. / = 1, ^ = 'uA'(n+ 

2 

(c?t, 9^, dyj, dx + vdw, 2tdt + xdx + vdy + 2wdw, fitdt + vdx-u{u + l)du-vdy + {^-2w)diy); 

9. f = l,A = u-^e^/": 

dx + vdw, 2tdt + xdx + vdv + 2wdw, tdt - vdx + v?du - ^dw). 

Equations with parameter-functions taking the values from cases 1 and 3 are locally equiv- 
alent to the constant coefficient equations from cases 2 and 4 correspondingly. Reducibility of 
equations with parameter- functions taking the values from cases 2, 4 and 5-9 and their potential 
symmetries to equations possessing Lie symmetries has been studied in [1,2,7]. 

This completes investigation of potential symmetries of diffusion-convection equations of 
form dS]) and consequently of form ([T]) . 

5 Exact solutions of Fujita— Storm equation 

One of the possible applications of the obtained potential symmetries and potential equivalence 
transformations is construction of exact solutions of equations ([T]) . A lot of examples of direct 
reductions of diffusion-convection equations with respect to potential symmetries can be found 
in [26-28]. Another way of finding exact solutions for equations of form ([T]) having nontrivial 
potential symmetries is to reconstruct them from the known solutions of potentially equivalent 
equations by means of application potential equivalence transformations. 
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Using the potential hodograph transformation (|16p we can reduce the Fujita-Storm equation 

Ut = {u~'^Ux)x (21) 

to the hnear heat equation v^: = Vxx ■ The same transformation connects exact solutions of these 
equations. 

Thus, e.g., the fundamental (source) solution v = (47rt)~"'^/^e~^ and dipole solution 
V = — ((47rt)~^/^e~^ ''^^*^)z of the linear heat equation are mapped into the separable and self- 
similar solutions of the Fujita-Storm equation 



u = (47rt)^/2gt.2 ^ ^ ^-1/2 / ^-y^cLy, 

Jo 

and „ = .-.(2f)'/^(l„^) 



correspondingly [3,21]. 

Applying the hodograph transformation to other known solutions of the linear heat equa- 
tion(collected, e.g., in [20]) we get the following exact solutions of the Fujita-Storm equation: 

u = c, u = x'\ u={x-2t)-'/\ u={x'±e^')-'/\ u = ±{e-'' - x')-'/^ 

1 t 
u = — ^ — , u 



4V24t2 + xv'-6t± V24t2 + x' x^-tln{xVi) 

1 



+ 2e-8* + 2e~^^^xJA - e8*(-cie-* ± + 2e-8* + 2e-'^^x) 



2 



Vcfe2* + 2e8* + 2e4<xVe-8*(-cie-* ± ^/cfe-^' + 2e-8* + 2e-4tx)2 - 4 



6 Conclusion 

This work is the fourth part of the series of papers on Lie group properties of nonlinear vari- 
able coefficient diffusion-convection equations of form ([T|). In the first part of the series [9] 
we construct the equivalence transformations and perform the complete group classification of 
class ([1]). We find not only the usual equivalence group G~ but also the extended one G"^ 
including transformations which are nonlocal with respect to arbitrary elements. Group G~ 
has a non-trivial subgroup of gauge equivalence transformations. In spite of really equivalence 
transformations, a role of gauges in group classification comes not to choice of representatives 
in equivalence classes but to choice of form of these representatives. Application of such gauge 
and nonlocal transformations is important for solving the problem of group classification in 
class (dJ. Moreover, let us note that for the class under consideration it seems impossible to 
obtain complete group classification in explicit form without using transformations from the 
extended equivalence group [9,12]. A similar statement is true for the problem of finding conser- 
vation laws: classification result with respect to the usual group of equivalence transformations 
can be formulated in an implicit form only. At the same time, using the extended equivalence 
group G"" , we can present the result of classification of conservation laws in a closed and simple 
form. 

The existence of the non-trivial subgroup of gauge equivalence transformations allows us 
to choose different gauges for simplification of the solving the group classification problem. 
The results of group classification for class ([1]) in two different gauges {g = 1 and g = h) are 
adduced in the first part [9]. In contrast to the most of papers on the subject under consideration 
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we performed two essentially different group classifications: classification with respect to the 
extended equivalence group and classification with respect to the set of all point transformations. 
The last classification shows the way to obtain exact solutions of quite complicated equations 
that are reducible to simpler form under action of point (additional) equivalence transformations. 
A number of examples of such solutions are adduced in the second part [10] of the series. We 
considered also in more details an interesting equation having sl{2, M)-symmetry and constructed 
a set of its exact solutions by means of Lie classical method, nonlinear separation of variables that 
is equivalent to a reduction with respect to a generalized nonclassical (conditional) symmetry 
generator. 

The next key point of the series is introducing the notion of contractions of equations, sym- 
metries, solutions, conservation laws and characteristics of conservation laws (the second [10] 
and the third parts [11] of the series). Investigation of contractions leads to deeper understand- 
ing connections between the different cases of classifications, it can be helpful for finding exact 
solutions and spaces of conservation laws of contracted equations directly from the ones of origin 
equations. 

The notion of equivalence of conservation laws with respect to a group of transformations has 
been generalized in several directions: classification of pairs "system -|- space of conservation 
laws"; classification of conservation laws for a given system with respect to its symmetry group; 
classification of pairs "system -|- a conservation law". Investigation of different generating sets 
of conservation laws of differential equations makes possible to generalize essentially procedure 
of construction of potential systems, and therefore, of finding potential symmetries. Namely, 
previously, for construction of simplest potential systems in cases when the dimension of the 
space of local conservation laws is greater than one, only basis conservation laws were used. 
However, the basis conservation laws may be equivalent with respect to groups of symmetry 
transformations, or vice versa, the number of C-independent linear combinations of conserva- 
tion laws may be greater than dimension of the space of conservation laws. The first possibility 
leads to an unnecessary, often cumbersome, investigation of equivalent systems, the second one 
makes possible missing a great number of inequivalent potential systems. Class ([I]) is very rich 
from this point of view. Namely, it contains nontrivial subclasses illustrating all these three 
possibilities. For all these subclasses we constructed the complete lists of inequivalent potential 
systems and investigated their potential symmetries. 
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